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Abstract. We show that coalescence of nucleons emitted prior to thermalization in highly excited nuclei can 
explain the anomaly of kinetic energies of helium clusters. A new coalescence algorithm has been included 
in the statistical approach to nuclear reactions formerly used to describe intermediate mass fragment 
production. 
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1 Introduction 



In this paper we address a phenomenon which seems to 
be important for understanding mechanisms which favour 
light cluster production in intermediate energy reactions . 
It concerns the so-called 3 He- 4 He "puzzle" , i.e. the anoma- 
lous behaviour of their kinetic energies. Usually, thermal 
models (e.g. [0) or approaches which consider a possi- 
ble radial expansion (flow) || are applied to describe the 
kinetic energies of light charged particles (LCP). In the 
thermal scenario we expect (Eki n (, 3 He)) (Ekini^He)) 
whereas radial flow delivers (Eki n ( 3 He)) < (Eki n { A He)). 
However, in many reactions quite an opposite behaviour 
has been observed. The corresponding data are summa- 
rized in table 1. For example, there is evidence of the "puz- 
zle" in the inclusive kinetic energy spectra of He isotopes 
obtained from p + Ag || and p+C [Q reactions at 1 GeV, 
in 7.5 GeV/c proton collisions with 12 C, 112 ' 124 Sn, 197 Au 
H and in Ne+U reactions at 250 and 400 A-MeV %. The 
anomaly was also observed in antiproton reactions 202 
MeV/c p + 12 C, 40 Ca, 63 Cu, 92 < 98 Mo and 238 U §. Precise 
measurements of 55 MeV 3 He+Ag collisions at <9; Q b=147.5° 
also have shown the anomalous behaviour of 3 He and 
4 He (ref.Q). The following trends have been observed in 
heavy-ion collisions: (i) the difference AE— (Ekmi 3 He))- 
(Ekin( 4 He)) is positive and increases with the particle 
multiplicity in Au + Au reactions at 250 A-MeV H, (ii) 
the " puzzle" is pronounced in central event samples of Au 
+ Au @ , 95 A-MeV Ar + Ni |n) , 50 A-MeV Xe+Sn @ 
but weaker in Ar + Ca Jll| collisions, (iii) positive values 
AE arc observed in the incident energy range from about 
50 A-MeV to 300 A-MeV but the anomaly vanishes at ~ 
1 GeV [Q. The existing data indicate that the larger the 
number of nucleons of the colliding system the larger the 
deviation of the kinetic energies of 3 He and 4 He. 



Table 1. Experimental data concerning the 3 He- 4 He anomaly. 
The temperatures T of the helium isotopes were obtained by 
Maxwell-Boltzmann fits to the kinetic energy distributions, in 
case of *) an exponential fit to the data at 6>; a b=90° was ap- 
plied. 



Reaction 


Incident energy 
or momentum 


T( 3 He) 
(MeV) 


T( 4 He) 
(MeV) 


p+ rz C\i\ 


1 GeV 


6.8+0.2 


5.8+0.2 


P+Ag g 




1 GeV 


11.2+0.2 


4.6+0.1 


P+ V2 C [1[ 




202 MeV/c 


19.8+12. 


15.6+9. 


P+ w Ca t 


1 


202 MeV/c 


24.2+13. 


16.8+9 


P+ M Cu [ 




202 MeV/c 


22.2+10. 


15.4+8. 


pV"Mo \\ 




202 MeV/c 


25.2+17. 


16.2+9. 


p+ MS Mo 




202 MeV/c 


21.8+6. 


17.1+7. 






202 MeV/c 


20.1+12. 


14.2+4. 


p+ ia C [5 




7.5 GeV/c 


37+2. 


28+2. 


p+'^Sn 5[ 


7.5 GeV/c 


38+1. 


27.il. 




7.5 GeV/c 


43+1. 


27.il. 


p+Au pi 


7.5 GeV/c 


50+2. 


33+2. 


2U Ne+U I 




250 A-MeV 


36.6+0.3* 


31.9+4.2* 


2U Ne+U 1 




400 A-MeV 


46.5+1.3* 


37.3+3.3* 






(MeV) 


(MeV) 


P+Ag H 


1 GeV 


17.7+0.3 


12.7+0.2 


^He+Agjs] 


55 MeV 


19.8 


15.6 


Ar+ Ni J 


J 


50-100 A -MeV 


32.5 


26.2 


Xe+Sn |1 


2] 


50 A-MeV 


51.0 


33.0 


Ar+Xe |1 




400 A-MeV 


139.6+2.5 


134.0+6.1 


Au+Au | 





100 A-MeV 


85.0+4.0 


59.0+2.0 


Au+Au [ 





150 A-MeV 


101.0+4.0 


77.0+5.0 


Au+Au [ 


of 


250 A-MeV 


147.0+4.0 


136.0+3.0 
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An appropriate way to describe processes involving 
many particles is the statistical approach. The system 
characterized in the initial stage by noncquilibrium dis- 
tribution functions evolves towards equilibration as a re- 
sult of interaction between particles. In this process the 
system runs through different states. The first one can 
be considered as equilibration of the one-particle degrees 
of freedom. In the following the evolution toward total 
thermalization can be considered as progressive involving 
of higher order particle correlations. For finite expanding 
systems one cannot predict in advance what kind of equi- 
libration should be considered. 

It is usually accepted that Intermediate Mass Frag- 
ments (IMF) observed in multifragmentation processes are 
mainly produced in a state close to thermalization. This 
assumption is supported by the success of statistical mul- 
tifragmentation models, e.g. SMM ||| and MMMC |||. 
However, this conclusion may not be true for composite 
particles like 3 He and 4 He which can be formed in earlier 
stages of the reaction. We guess that their unexpected be- 
haviour results from an interplay of different production 
mechanisms. 



2 Production of composite particles by 
coalescence 

First we point at an alternative way to form clusters. We 
start from a distribution of nucleons in the phase space at 
some "freeze-out" time obtained after a dynamical evolu- 
tion. Generally, non-uniform distributions are conceivable. 
But in some experiments, e.g. central nucleus-nucleus col- 
lisions jl7j , it is possible to select events which are nearly 
isotropic in space and look like thermal ones. Therefore, in 
such cases, we can simply assume that the nucleons pop- 
ulate the available many-body phase space uniformly, i.e. 
there is equilibration in one-particle degrees of freedom, 
that gives rise to a thermal distribution for individual nu- 
cleons in thermodynamical limit. 

A composite particle can be formed from two or more 
nucleons if they are close to each other in the phase space. 
This simple prescription is known as coalescence model 
and it reflects the properties of the nucleon-nucleon inter- 
action. Here we use the coalescence in momentum space 
only. The basic assumption is that the dynamical process 
which leads to a momentum redistribution is very fast 
(nearly instantaneous) so that the coordinates of nucleons 
are just defined by their momenta. It is also justified tak- 
ing into account quantum properties of the system since 
the wave functions of nucleons are rather broad. This type 
of coalescence model has proven successful in reproducing 
experimental data (see e.g. |ig|,|i1|p0|]). 

In the standard formulation of the model it is assumed 
that the fragment density in momentum space is propor- 
tional to the nucleon density times the probability of find- 
ing nucleons within a small sphere of the coalescence ra- 
dius pq. From this hypothesis an analytical expression can 
be derived for differential yields of coalescent clusters 



@0: 

V dp\ ) 2 A v\z\\i m n J V d Pn J 

where Ea, E n ,pA and p n are the energies and the mo- 
menta of fragments and nucleons, v and z are the numbers 
of neutrons and protons [A = i/ + z); Sa is the spin of the 
fragment with mass number A and m n m 0.94 GeV is the 
nucleon mass. However, this equation disregards correla- 
tions between different clusters since the conservation of 
the nucleon number is not taken into account. Therefore, 
the above formulae is valid only for Na=i 3> Na=2 S> 
Na=3--- In many reactions this condition is not fulfilled. 

We developed an alternative formulation of the coales- 
cence model which is suitable for computer simulations. 
Nucleons can produce a cluster with mass number A if 
their momenta relative to the center-of-mass moment of 
the cluster is less thanp - Accordingly we take |pj— p cm | < 

Po for all i = I,..., A, where p cm = -jYst^iP*- Tnis is 
performed by comparing the momenta of all nucleons. In 
the following examples the value po « 94MeV/c has been 
adopted corresponding to relative velocities v re i=0.1c in 
agreement with previous analyses |jl9|| . 

We note a problem which is sometimes disregarded in 
these simulations. Some nucleons may have such momenta 
that they can belong to different coalescent clusters ac- 
cording to the coalescence criterion. In these cases the final 
decision depends on the sequence of nucleons within the 
algorithm. To avoid this uncertainty we developed an iter- 
ative coalescence procedure. M steps are calculated in the 
coalescence routine with the radius poj which is increased 
at each step j: poj = (j/M) ■ po (j = 1,...,M). Clusters 
produced at earlier steps participate as a whole in the fol- 
lowing steps. In this case the final clusters not only meet 
the coalescence criterion but also the nucleons have the 
minimum distance in the momentum space. Mathemati- 
cally exact, this procedure gives correct results in the limit 
M — > oo but we found that in practical calculations it is 
sufficient to confine the steps to M= 5. 

To demonstrate how fragments are produced by co- 
alescence we take as an example a nuclear system with 
mass number Ao = 200 and charge Zq = 80 at various ex- 
citation energies E*=1Q, 20, 30 and 100 A-MeV. We dis- 
integrate the system into nucleons by taking away about 
7 A-MeV (binding energy). The rest of the energy turns 
into the kinetic energies of nucleons which populate the 
whole available many-body momentum phase space uni- 
formly. We use the procedure developed in pl[ to generate 
momenta. Example mass distributions of generated frag- 
ments are shown in fig. [lj At low excitation energies we see 
one big " coalescent" cluster and some light fragments (an 
"U-shape" distribution). With increasing excitation en- 
ergy the large cluster will be destroyed and several small 
clusters are produced. At about 30 A-MeV a transition ap- 
pears which changes the mass distribution into a "step"- 
like one. At higher energies an exponential-like shape of 
the mass distribution occurs and that corresponds to the 
standard coalescence picture at high energies. We found 
that this transition happens at (p)/po ~ 2, where (p) is 
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Fig. 1. Fragment production by coalescence. Upper panel: 
fragment yields per 10 4 generated events at 4 different exci- 
tation energies. Lower panel: the corresponding fragment mo- 
menta per nucleon. 



the average momentum of thermalized nucleons, and it 
depends weakly on the size and other parameters of the 
system. In fig. [j] we show also the average momenta of 
the clusters. We see that the coalescence provides a spe- 
cial ordering in the kinetic energy of the fragments: light 
clusters have larger energies per nucleon than the heavy 
ones. Large clusters are preferentially formed from nucle- 
ons with small momenta, i.e. which are closer to the center 
of the phase space. 

The evolution of the mass distribution with excitation 
energy predicted by the coalescence model corresponds to 
our expectations for the decay of finite nuclei. Large clus- 
ters must be excited if they are treated in the classical coa- 
lescence model because there is a motion of nucleons with 
respect to the center of mass of the clusters. In this case 
a big coalescent cluster can be treated as an excited rem- 
nant of the system in which more complete equilibration is 
achieved and which can decay in a thermal-like way. Such a 
picture is supported by dynamical calculations which show 
the different reaction stages in nucleus nucleus collisions 
at higher energies (e.g. P2fl ). After the first interactions 
some nucleons immediately leave the system and we can 
assume that later on they form coalescent clusters in final 
state interactions. This fast emission lowers essentially the 
energy of the remaining part and favours its equilibration, 
so we have a chance to treat it by statistical methods. 



1800 



1500 



thermal 



E =20A.MeV 
A =200, Z =80 




coalescent 



100 



> 75 

A 



150 200 

E kin (MeV) 



50 



LU 

V 



; v 


200, Z =80 
















coalescent 








thermal - 











25 



50 



75 

E (A.MeV) 



100 



Fig. 2. Kinetic energy spectra (top) and mean kinetic energies 
per nucleon (versus excitation energy, bottom) of 3 He and 4 He 
calculated for the thermal and coalescence scenaria. 



3 Kinetic energies of 3 He and 4 He 

Now we show qualitatively how the coalescence can ex- 
plain the " puzzle" . In fig. ^ (top) we present the kinetic 
energy distributions of coalescent clusters 3 He and 4 He 
generated for E* = 20A-MeV. The peaks in both spec- 
tra are shifted towards higher energies similar to collective 
motion (radial flow). For a different (e.g. non-equilibrium) 
initial distribution of the nucleons in momentum space we 
expect broader spectra with possible shifts of the peaks. 
But, in any case, the kinetic energies of the clusters will be 
high. In the bottom of fig. 2 we plot the mean kinetic en- 
ergies of helium isotopes for two scenaria and different ex- 
citation energies. The coalescent 3 He clusters have slightly 
larger kinetic energy per nucleon compared with 4 He. This 
energy difference increases with excitation energy though 
the total kinetic energy of 3 He is still lower than the en- 
ergy of 4 He. This is caused by the mentioned correlation 
between the cluster production and the momenta of nu- 
cleons which form the clusters (fig. 1, bottom). Such a 
correlation was already under discussion, e.g. y23|, to ex- 
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plain the observed radial-flow energies of the fragments. 
In fig. 2 we also show the kinetic energies and spectra of 
3 He and 4 He calculated with the SMM code assuming 
a completely thermalized system at the same excitation 
energies. Herewith all thermal processes including deex- 
citation of hot primary fragments are taken into account. 
There is a striking difference in the kinetics of light clusters 
produced by coalescence or thermal mechanisms. In high 
energy reactions the difference may even increase because 
the energy released in the preequilibrium stage is usually 
much larger than the remaining equilibrium energy. We 
have checked also the influence of secondary deexcitations 
of possibly excited intermediate coalescent clusters (with 
A > 4) on the spectra of the helium isotopes: the energy 
decrease is relatively small and the difference to the ther- 
mal spectra remains pronounced. 

We expect that the thermal fragment production and 
the coalescence coexists in heavy-ion collisions over a wide 
energy range. The final shapes of the fragment spectra 
depend on their relative contributions. It is well known 
that the equilibrium decay provides a large number of a- 
particles, whereas 3 He clusters are suppressed. But 3 He 
clusters can be produced predominantly by coalescence 
from the large number of emitted preequilibrium nucle- 
ons. As a result, the kinetic energies of 3 He become larger 
compared with those of the ev-clusters. Such conditions 
may occur at intermediate energies in heavy ion reactions. 
However, in central nucleus- nucleus collisions the size of 
the thermalized source decreases with increasing beam en- 
ergy as shown by SMM calculations and comparison with 
experimental data |^,^5). Thus at sufficient high energies 
the coalescence becomes a dominant production mecha- 
nism for both kinds of He clusters and their kinetic ener- 
gies may be very similar. 
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Fig. 3. Mean kinetic energies vs. mass number of LCP's in 
central collisions of Au on Au at incident energies 250 and 
1150 A-MeV. Black symbols: data, dotted line: thermal pro- 
duction only, shadowed area: thermal orign and coalescence 
(v re i=0.1c) assuming an overall uncertainty of ~ 15 % in the 
model parameters. 



4 Comparison of experimental data with 
SMM and coalescence calculations 

We checked this hypothesis by calculations of fragment 
production in central Au + Au collisions using the SMM 
model which was modified by inclusion of radial flow and 
preequilibrium particle emission. Recently this approach 
was used in ref's.|2^,|4],[2^,|2(|. In the framework of this 
model it is assumed that the total cm. energy is shared be- 
tween a thermal source (responsible for IMF and thermal 
LPC production), creation of pions (taken phenomenolog- 
ically from experimental data) and preequilibrium emis- 
sion of nucleons. Afterwards, these nucleons may coalesce 
into composite particles. A uniform initial distribution 
of "preequilibrium" nucleons in many-body momentum 
space is assumed. As already mentioned this assumption 
allows to estimate simply the influence of coalescence in 
the particular case of central events. 

The excitation energy and mass of the equilibrium 
sources were parametrized by fitting the charge distri- 
butions and multiplicities of IMF's in the beam energy 
range from 150 to 1050 AMeV as described in ref. fl24f . 
The radial flow was found by fitting the kinetic energy 



spectra of IMF's and it matches the experimental sys- 
tematics |l7|,[27j. The source parameters at 1150 AMeV 
[[Til obtained by extrapolation are very close to that at 
1050 A-MeV. In ref's. (HH) it has been concluded that 
the size of the thermal source decreases considerably with 
the beam energy, and preequilibrium emission becomes 
important. The part of energy attributed to the preequi- 
librium stage is very large and, therefore, the masses of 
the coalescence fragments decrease nearly exponentially. 
In fig. ^ we compare the experimental data from central 
Au+Au collisions with the calculated mean kinetic ener- 
gies of light charged particles. The energy spectra of 3 He 
and 4 He were simulated for 250 AMeV in the angular 
range 60° < <9 cm < 90° using the lower and upper reg- 
istration thresholds given in ref.jl0|]. LCP's emitted from 
the thermal source with superimposed collective motion 
show the expected increase of (Eki n ) with A. However, 
the energy stored in preequilibrium nucleons is decisive 
for clusters predominantly produced from these nucleons: 
the 3 He- 4 He anomaly appears if a substantial part of 3 He 
comes from coalescence whereas the thermal production 
of 4 He dominates. In addition we show in table 2 a rea- 
sonable reproduction of the hydrogen and helium multi- 
plicities within this scenario. 
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Table 2. Multiplicities of LCP's in central Au+Au collisions. 
Typical sum uncertainties of the calculations are ± 5. 
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Fig. 4. Yield ratio 3 He/ 4 He in central Au+Au collisions vs. 
beam energy. Symbols show the data, dotted line: equilibrium- 
source calculations with the parameter set from p4[], shad- 
owed area: calculations including preequilibrium and coales- 
cence, similar to fig. 3. 



Enhanced production of 3 He with increasing beam en- 

ty has been established in central Au+Au collisions jlo|, 
p9[ as shown in fig. |j. The 3 He yield dominates over 
ne production in the region > 600 A-MeV. Calculations 
with our model confirm this trend. There is an obvious 
surplus of 3 He by coalescence with respect to the thermal 
production. This behaviour seems to be related to the dis- 
appearance of the anomaly established at 1.15 A-GeV p4[ , 
fig. ||, that is also present in our calculations. The disagree- 
ment in the calculated proton energies may be caused by 
the centralitycriterion in ref. E3 which is different from 



that in ref. 10 . Selection of protons around cm ~ 90° 
in events with high multiplicity may push down their ki- 
netic energies since such a criterion includes many events 
with predominantly forward-backward directed high en- 
ergy nucleons. Indeed, lower proton energies were obtained 
by simulations assuming anisotropic preequilibrium emis- 
sion of nucleons whereas the conclusion about the 3 He- 
4 He-anomaly is maintained. 



5 Conclusions 

Summarizing we found that the appearance and disap- 
pearance of the 3 He- 4 He anomaly reflects an interplay 
of equilibrium and nonequilibrium processes in these re- 
actions. Besides the thermal emission the coalescence of 
nonequilibrium nucleons was found to be the important 
process which is responsible for enhanced mean kinetic 
energies of 3 He. The effect is pronounced when both the 
thermal production and the coalescence contributions are 
comparable. At sufficient high energies only a small part 





E(A-MeV) 


pre 


equ 


sum 


data, ref. jl7| 




150 


32.2 


32.4 


64.6 


61.84(0.58) 


Z=l 


250 


52.2 


30.8 


83.0 


75.82(0.62) 


Z=l 


400 


60.0 


34.8 


94.8 


92.04(0.62) 


Z=2 


150 


5.3 


15.6 


20.9 


26.76(0.36) 


Z=2 


250 


8.9 


14.5 


23.4 


27.27(0.36) 


Z=2 


400 


9.3 


14.3 


23.6 


24.16(0.30) 



of the colliding nuclei attains equilibration and both 4 He 
and 3 He are mainly produced by coalescence. As a conse- 
quence, the anomalous behaviour vanishes. 
We mention another approach ||l2f which shows that in 
central nucleus- nucleus collisions around the Fermi energy 
a fast fragment emission during the expansion stage (de- 
scribed by the EES model (3^] ) can be responsible for the 
" puzzle" . This explanation is qualitatively in agreement 
with our findings. However, we suppose that coalescence 
is a general mechanism which may be applied successfully 
in statistical or dynamical models (see e.g. |^,^l],|3^]). It 
seems that the correlation between energies and masses 
of the clusters predicted by the coalescence model may 
explain the observed features of fragment production in a 
wide energy range. 
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